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We apply ab-initio quantum chemical methods to calcu- 
late correlation effects on cohesive properties of NiO, thereby 
extending a recently proposed scheme to transition metal ox- 
ides with partially filled d-bands. We obtain good agreement 
with experiment for the cohesive energy and show that the 
deviation of the lattice constant at the Hartree-Fock level is 
mainly due to van der Waals-like interactions. Correlations 
enhance the stability of the magnetic ground state found at 
the Hartree-Fock level. 



containing transition-metal compounds. In addition, we 
want to compare some relevant atomic and molecular re- 
sults from ab-initio quantum chemistry to those of DFT. 
For this purpose, we performed test calculations with the 
local density approximation (LDA) as-well as with the 
gradient-corrected functional of BeckeEll, Lee, Yang and 
Parr (BLYP)E3. 



II. THE METHOD 



A. Incremental scheme 
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I. INTRODUCTION 

NiO has a long tradition in solid-state physics as-a pro- 
totype of an insulator with partially filled rf-bandad. The 
initial, explanation for its unusual behavior was given by 
MottH who argued that the Coulomb repulsion of the elec- 
trons in singly occupied nickel orbitals should completely 
suppress excitations of the d — > d type. This picture led 
to the notion ofJVlott insulators and was later refined 
by Zaanen et alQ who characterized NiO as a charge- 
transfer insulator with a valence band of oxygen p-type 
and a resulting p — ► d gap. The band structure and es- 
pecially the gap have been in the center of interest and 
a considerable-appunt of work has been devoted to their 
determinationaEJ (for reviews on NiO see, e.g., Refs. 
[l5|,|l6] and, more general, Ref. Surprisingly, only lit- 
tle work has been done on cohesive properties and we are 
aware pi, only one result from density functional theory 
(DFT)lia. At the Hartree-Fock (HF) level, a very-care- 
ful and extensive investigation has been performeco. In 
the present paper, we want to address the question how 
electron correlations affect ground-state properties. 

In the field of atoms and molecules, the most accu- 
rate methods are those of quantum chemistry whereas in 
solid-state physics the dominant method is DFT. How- 
ever, it has been demonstrated by using local opera- 
tors combined with an incremental method that high- 
quality results can be obtained, from a quantum-chemipa.1 
approach to semiconductorsEfl and ionic compoundsEl 
Transition metals are still a challenge in quantum chem- 
istry and even for the NiO molecule extensive calcula- 
tions are necessary to obtain agreement with experiment. 
We want to show, however, that an accurate determi- 
nation of correlation effects is possible nowadays, using 
the afore-mentioned incremental scheme, even for solids 



The incremental scheme uses the fact that electron cor- 
relations are a local property. The first step is to calcu- 
late the correlation energy e(A) of groups -of localized 
orbitals (the so-called one-body increments)Ei The cor- 
relation energy of two groups of orbitals e(AB) can then 
be expressed as 

e(AB) = e(A) + e(B) + Ae(AB) , 

the non-additive part Ae(AB) being the two-body incre- 
ment. For three groups of orbitals one obtains 

e(ABC) = e(A) + e(B) + e(C) + 
Ae{AB) + Ae(BC) + Ae{CA) + Ae(ABC). 

Obviously, this series can be extended to any order of 
increments and for the correlation energy of the solid 
one obtains 



Cbulk 



E e ^) + ?E Ae (^) + ^ E m^bc) 



A.B 



A,B,C 



(1) 



which is still exact. The approximation consists in ne- 
glecting higher than third orders of increments which 
makes it necessary that the sum of these increments 
should be negligibly small. Moreover, it is assumed 
that the increments rapidly decrease with increasing dis- 
tance between groups A and B. Since the calculations are 
performed on embedded clusters, the increments should 
only weakly depend on the chosen cluster. The valid- 
ity of these assumptions has to be checked, of course, 
for each specific case. A prerequisite is that the corre- 
lation method must be size-extensive. Since the method 
must be applicable to all types of wavefunctions (includ- 
ing low-spin), we chose the quasidegenerate variational 
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perturbation theory In the simplest case of 

a single-reference wavefunction, the correlation energy is 
obtained from the functional 



, = < ^scf + ^c\H - E scf \^scf + 9 
scf\^scf > 



> 



where ^scf is the self-consistent field wavefunction and 
"fc the correlation part obtained by single and double 
substitutions from the SCF determinant, i.e., 



I*c >= 



z r a a+a a \$>sCF > + c 7 a s b a+ a+ a a a b \^> scf > 



The correlated orbitals are 2s and 2p for O, 2 ~ and 
the 3c?-orbitals in the case of Ni 2+ . The doubly occupied 
orbitals were kept closed in the reference wavefunction 
which leads to a single reference calculation except for 
the case of low-spin coupling. The cohesive energy was 
calculated with respect to the Ni 3d 8 4s 2 state (of course, 
the 4s orbital was correlated in that case, too). The in- 
fluence of the correlation of the 3s and ip orbitals of Ni 
will be discussed separately. 



C. Test calculations 



The method is closely related to the coupled electron pair 
approximation (CEPA(Ql)cJ and the averaged coupled- 
pair functional (ACPF)Ej, see the discussion in Ref. |27]. 
The largest increments were compared with those ob- 
tained from the coupled-cluster .approach with single 
and double substitutions (CCSD)Ea where an exponen- 
tial ansatz for the correlated wavefunction is made: 

\^ccsd >= exp(^ c r a a+a a + ^ c r a s b a^afaaab)\9scF > 



These increments were also calculated with CCSD(T)@ 
where three particle excitations are included by means of 
perturbation theory. All the calculations [ Were d one with 
the ab-initio program package MOLPROlEIlIIi 2 -! 

B. Pseudopotentials and basis sets 

The calculations are performed on clusters of two or 
three ions which are accurately treated by quantum- 
chemical methods, embedded in surroundings treated at 
a lower level of approximation. To model the Pauli re- 
pulsion of the surroundings on the clusters considered, 
a pseudopotential for Ni 2+ was employed simulating a 
fixed [Ar] 3d 8 configuration, cf. Tabic |. The parame- 
ters have been adjusted to averaged valence energies of 
12 configurations of Ni and Ni + . The valence energies 
were taken from relativistic HF calculations using the 
Wood-Boring HamiltonianE 2 ]. The clusters were modeled 
in the same way as described in previous work on MgO 
and CaOEL Cluster O 2 " ions were surrounded by Ni 2+ 
pseudopotentials as next neighbors. The whole system 
was embedded in a set of point charges (typically 7 x 
7x7 lattice sites with charges ±2 in the interior and 
reduced by factors of 2, 4 and 8 at the surface planes, 
edges and corners, respectively). The description of the 
ions contained within the clusters is aa_follows. The oxy- 
gen basis set was a [5s4p3d2/]-basista. For the nickel 
ions a pseudopotential with a Ne coreE-3 together with an 
optimized basis set described in Ref. |34| was used. The 
d-contraction was slightly relaxed (the three smallest ex- 
ponents were used uncontracted) and three /-exponents 
(0.67,2.4,7.8) were added yielding a [6s5p4d3/] basis set. 



1. Nickel ionization potential 

Results of test calculations for the first and second ion- 
ization potential of Ni are shown in Table O. We per- 
formed calculations with a finite basis set (indicated as 
'our basis set') as well as, when possible, using numerical 
schemes which are comparable to a complete basis set (in- 
dicated as 'finite differences'). Comparing Hartree-Fock 
and Dirac-Fock calculations, one finds that relativistic 
effects make a significant difference. The origin and im- 
portance of these effects on the results for the crystal 
binding energy will be discussed later. The HF calcula- 
tion with the Ne-core pseudopotential and the [6s5p4d3/] 
basis set is in nice agreement with the Dirac-Fock result 
which shows that the pseudopotential is well suited to ac- 
count for relativistic corrections. At the correlated level, 
we find that we recover most of the missing energy. Cor- 
relation of Ni 3s and 3p orbitals slightly improves the 
result; the remaining error is mostly due to basis-set de- 
ficiencies. The density functional results tend towards an 
overestimation. 



2. Electron affinity of oxygen 

The next test concerns the electron affinity of the oxy- 
gen atom. We find a-j/alue of 1.359 eV at the QDVPT 
level (exp.: 1.461 eV£3). Note that with several density 
functionals negatively charged ions are not stable because 
of the wrong asymptotic behavior of the functionals (see 
the discussions in Ref. |36|-p8|) . Indeed, we found positive 
orbital energies when applying LDA and BLYP to the 
0~ ion which indicates the instability. 



3. Ground state of the NiO molecule 

A final test calculation was performed for the 3 S _ 
ground state of the NiO molecule (see Table |||). Our 
multi-reference averaged coupled-pair functional (MR- 
ACPF) result for D e is in good agreement with the ex- 
perimental result. The bond length is slightly too short 
and therefore the vibrational frequency is too high (cal- 
culating the vibrational frequency at the experimental 
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bond length gives a value of 897 cm" 1 which is in better 
agreement with experiment) . Enlarging the basis sets by 
using g functions leads to a further slight reduction of the 
bond length. Correlating the Ni 2>p orbitals also reduces 
the bond length (by ~ 0.005 A). Therefore, an enlarge- 
ment of the active space (the group of orbitals from which 
the reference determinants of the MR-ACPF calculation 
are formed) seems to be necessary. The lowest-lying vir- 
tual orbitals not yet included in the active space are of tt x 
and TT y type. However, inclusion of these orbitals is not 
feasible currently due to the high computational effort. 
Bauschlicherc3 added another virtual a orbital to the ac- 
tive space, but on the other hand the 5 orbitals and the 
O 2s-orbital were kept closed (i.e. correlated, but dou- 
bly occupied in the reference wavefunction) . This gave 
excellent agreement with experiment. However, it is ob- 
vious from his work that the results are still sensitive to 
the correlation level applied (see the comparison in Ref. 
|39|) . Actually, we found that this problem is more crit- 
ical in the case of the molecule than in the calculations 
performed for the solid (QDVPT calculations, e.g., are 
difficult to converge for the molecule). Note that, in con- 
trast to the solid, the NiO molecule is covalently bonded 
with charges ± 0.7 according to a Mulliken population 
analysis. The gradient-corrected DFT results for NiO 
are excellent for the vibrational frequency and the bond 
length, but for the dissociation energy one finds a sig- 
nificant overbinding. This might change with improved 
functionals such as B3LYP or B3PW91 which are not yet 
implemented in MOLPRO. 



III. RESULTS FOR THE SOLID 

A. Incremental expansion and magnetic properties 

Regarding the one-body correlation energy increments 
(see Table IV), one finds nearly the same value for the 
increment O — > 2 ~ as in the systems MgO and CaO 
whereas the increment Ni — > Ni 2+ is significantly larger 
than the coEtesponding value in the case of the alkaline 
earth oxidesc3. As compared to Mg, the low lying d- 
orbitals led to an enhanced core-valence correlation and 
a larger increment for Ca — > Ca 2+ . The near-degeneracy 
effects in the Ni atom lead to an even higher valence 
correlation energy. The two-body increments between 
two oxygen ions are nearly 50 % larger than those of MgO 
although the lattice constants are very similar. This may 



be connected with an increased polarizability of O in 
NiO which in turn might be related to the availability of 
low- lying open-shell d-orbitals in Ni 2+ . 

As in the systems MgO and CaO, the Mulliken popu- 
lation analysis gives anX^O charge transfer of slightly 
less than 2. The Ni-0 two-body increments turn out to 
make the most significant correlation contribution. This 
is not surprising since the polarizability of the open-shell 
Ni 2+ ion is roughly six times higher than that of Mg 2+ . 



We used a special procedure to evaluate the Ni-0 in- 
crements for technical reasons. Single external correla- 
tion contributions to e(0 2_ ) involving single substitu- 
tions from localized 2 ~ 2s, 2p orbitals into unoccupied 
orbitals, were not directly calculated for the embedded 
Ni 2+ 2 ~ cluster but were transferred from a calculation 
for a single O 2- embedded in large-core Ni 2+ pseudopo- 
tentials. 

Analyzing the correlated wavefunction in the case of 
the Ni-0 increments for next neighbors, one finds that in- 
ternal excitations (excitations from the oxygen ion to the 
singly occupied nickel orbitals) are negligibly small and 
contribute with a few /iH (Hartree units, 1 H=27.2114 
eV) only. This is certainly a consequence of the large 
Coulomb repulsion U in the case of doubly occupied 
Ni 3c?-orbitals. Another interesting feature is the im- 
portance of spin-flips on the Ni site. By restricting the 
excitations for the supersystem, the Ni 2+ and O 2- ions 
to double externals (only two-particle excitations to com- 
pletely unoccupied orbitals are allowed, excitations to or- 
bitals already occupied with an electron of opposite spin 
are forbidden), one obtains an increment of only 0.008450 
x6 H instead of 0.011534 x 6 H when single external 
excitations are additionally taken into account. Thus, 
the importance of single externals becomes obvious. The 
processes among the single externals which are the most 
important ones for the increment are those where an elec- 
tron from the O 2- ion is excited to a singly occupied Ni 
orbital and the Ni electron is excited to a virtual orbital 
which results in an effective spin-flip (Fig. [j]). Of course 
these processes are only possible in open shell systems. In 
the case of the oxygen-oxygen increment for next neigh- 
bors, neglect of single externals only slightly reduces the 
increment from 0.003916 x 6 H to 0.003549 x 6 H. This 
confirms that single externals have their main impact on 
increments involving open shells. 

The Ni 2+ -Ni 2+ increments were calculated for the case 
of high (S=2) and low spin (5*=0) coupling (neglecting 
single externals). However, the very small energy split- 
ting between the states is sensitive to the chosen cluster 
and increases when point charges as next neighbors are 
replaced by explicitly treated oxygens. Therefore, the 
results in Table IV should not be used for a quantita- 
tive evaluation of the splitting. In the case of a cluster 
of the type (O(0,0,0)-Ni(0,0,l)-Ni(0,0,-1)), the splitting 
is 138 /zH at the HF level and 276 /iH at the correlated 
level including single external excitations (the S=0 state 
is lower). For next neighbors, one finds a tendency to- 
wards ferromagnetic coupling: in the cluster (0(0,0,0)- 
Ni(0,0,l)-Ni(0,l,0)), the S=2 state is 19 lower at 
the HF level and 22 /iH when correlations are included. 
These results are in qualitative agreement with experi- 
ment where the exchange coupling is ferromagnetic for 
next neighbors and antiferromagnetic for second next 
neighborsEI This indicates that the stability of the an- 
tiferromagnetic AF2 state (a state with parallel spins in 
(1 1 1) planes, adjacent planes having antiparallel spin) 
found at the HF level, will be enhanced due to corre- 



3 



lations. Although quantitative values for the exchange 
coupling could be calculated for molecules and are in 
agreement with experiment (see, e.g., Ref. 41), the situa- 
tion is more difficult in solids and the exchange cou pl i ng 
increases with the number of ions explicitly treatedcHlc3. 
Therefore, our results can only give a qualitative picture 
here. Note that the order of magnitude of the magnetic 
splitting (0.03 eV/ double cell at the HF level) is much 
smaller than the total cohesive energy. 

The incremental expansion for the correlation contri- 
bution to the cohesive energy is well convergent. The 
two-body increments show a rapid van der Waals like 
decrease and the sum of the three-body corrections is 
small compared to the sum of one- and two-body incre- 
ments. It is interesting to see that the two-body incre- 
ments amount to roughly 80% of the correlation contri- 
bution to the cohesive energy, while the intra-ionic contri- 
butions contained in the one-body increments yield only 
~ 20% (in the case of the alkaline earth systems, both 
contributions were roughly of the same size). The QD- 
VPT results are closer to CCSD(T) than to CCSD, but 
the results are weakly dependent only on the method ap- 
plied. Individual increments and their total contribution 
to the bulk cohesive energy for a lattice spacing of 4.17 
A are displayed in Table IV and Figure 0. 



B. Comparison with experiment 

Our reference value for the HF binding energy (6.2 eV) 
is taken from Ref. 
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However, two corrections are nec- 
essary: since the calculations of Ref. |l4| were done within 
the unrestricted Hartree-Fock (UHF) approximation for 
the solid, but at the restricted Hartree-Fock (RHF) level 
for the atomsC-j, a part of the correlation energy is al- 
ready included in that value. We corrected this by sub- 
tracting the energy difference \Ejjhf — Erhf\ for the 
Ni 2+ ion (-0.0072 H) from the binding energy of Ref. 
[Ti] . This correction can be viewed as a 'one-body cor- 
rection' which is tantamount to applying Eq. [I] not to 
the correlation energy but to the difference UHF- RHF. 
Since the correction is nearly the same for the case of a 
system with a Ni 2+ and an 2 ~ ion, we conclude that 
the 'two-body corrections' are negligible and the HF co- 
hesive energy of Ref. H should be reduced by 0.007 H. 
A second correction is necessary, since the CRYSTAL 
calculation was done neglecting effects arising from rela- 
tivity. These effects gain importance beginning with the 
first row of transition metalscS. We found (see Table ||) 
that the first and second ionization potentials of Ni in- 
crease from 0.840 H (Hartree-Fock calculation) to 0.853 
H (Dirac-Fock calculation) . This is a consequence of the 
relativistic contraction (and stabilization) of the Ni 4s- 
orbital. In the case of Ca, this correction is of the order 
of 2.4 mH only. Spin-orbit effects are negligible for Ni 
as the d-occupation in the solid is d 8 (t^gSg) and very 
similar to that of the isolated atom ( 3 F). Taking both 



corrections into account, we obtain a cohesive energy of 
5.65 eV (5.59 eV at 4.17 A, using the bulk modulus from 
Ref. [l4j). Adding the zero-point vibrations (0.12 eV, a 
Debye approximation with 9 ~ 640-KL3 was used) to the 
experimentally determined cohesive energy of 9.5 eVE3, 
we thus deduce an 'experimental' correlation contribu- 
tion to the cohesive energy of « 4.0 eV. Within the basis 
sets used, we recover 3.38 eV (= 84 %) of this energy 
with our local correlation scheme. 



C. Lattice constant 

By performing the incremental expansion for a second 
value of the lattice constant a, we obtain the correla- 
tion correction to a. The same physical picture as in 
the case of MgO and CaO is found: the van der Waals 
like interaction between the ions (included in the two- 
body increments) reduces the lattice constant whereas 
correlations in the oxygen ion tend to enhance the lattice 
constant. The latter effect is due to the lower level spac- 
ing and therefore increasing importance of correlations 
at larger lattice constant. The results for a=4.264 A are 
displayed in Table |v|. Concerning the three-body incre- 
ments, we used the two largest ones to evaluate the lattice 
constant. By linearly fitting the correlation energy aud, 
using the bulk modulus from the CRYSTAL calculationEa 
(214 GPa), we obtain a value of 4.167 A, i.e., a reduction 
by 0.097 A with respect to the CRYSTAL-HF result. 
This is in-aiice agreement with the experimental result 
of 4.17 Acll. We assumed a NaCl structure and did not 
take into account the slight rhombohedral distortion ob- 
served in the experiment. The calculated cohesive energy 
at this equilibrium constant is (5.65-0.06+3.38) eV=8.97 
eV which is 93 % of the total cohesive energy of 9.6 eV. 

Finally, we evaluated the one-body and two-body Ni-0 
increments when electrons in the Ni 3s and 3p orbitals 
are correlated. At a lattice constant of 4.17 A, we find 
that the one-body increment Ni — > Ni 2+ increases from 
0.074545 to 0.080790 as a consequence of the outer-core- 
valence correlation. On the other hand, the absolute 
value of the sum of the Ni-0 increments increases from 
0.073446 to 0.082706. As these effects have the opposite 
sign, they nearly cancel and increase the correlation con- 
tribution to the cohesive energy to 0.127204 Hat 4.17 A. 
The lattice constant is reduced to 4.151 A. 



D. Comparison with results from literature 

Our results for the cohesive energy are in better agree- 
ment with experiment than those of DFT, where a signif- 
icant overestimation was foundlia. The lattice constant of 
the latter work (4.183 A) is in good agreement with ex- 
periment. A recently performed Quantum Monte-Carlo 
(QMC) calculatiorO gives good agreement with exper- 
iment, but the error-bars are quite large compared to 
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those of other QMC worked. Moreover, the basis sets are 
very small and not adjusted to the pseudopotential used. 



IV. CONCLUSION 

In conclusion, we have shown that the incremental 
scheme is capable of accurately treating systems with 
open d-shells. We obtained more than 84 % of the corre- 
lation contribution to the cohesive energy or 93 % of the 
total cohesive energy. Contributions from interatomic 
origin are roughly four times larger than those due to 
the formation of ions. The lattice constant agrees well 
with experiment. The stability of the magnetic ground 
state obtained found at the HF level is enhanced with 
the inclusion of correlations. Wc feel that the quality of 
this calculation is superior to previous work with DFT 
and QMC. 



ACKNOWLEDGMENTS 

We are grateful to Prof. H.-J. Werner (Stuttgart) for 
providing the program package MOLPRO and to Dr. N. 
M. Harrison (Daresbury) for useful information on the 
CRYSTAL calculations. 



1 J. H. de Boer and E. J. W. Verwey, Proc. Phys. Soc. 49 
(extra), 59 (1937) 

2 N. F. Mott, Proc. Phys. Soc. London A 62, 416 (1949) 

3 J. Zaanen, G. A. Sawatzky, J. W. Allen, Phys. Rev. Lett. 
55, 418 (1985) 

4 L. F. Mattheiss, Phys. Rev. B 5, 280 (1972); Phys. Rev. B 
5, 306 (1972) 

5 K. Terakura, A. R. Williams, T. Oguchi, J. Kiibler, Phys. 
Rev. Lett. 52, 1830 (1984) 

6 A. Fujimori, F. Minami, Phys. Rev. B 30, 957 (1984) 

7 G. J. M. Janssen, W. C. Nieuwpoort, Phys. Rev. B 38, 
3449 (1988) 

8 A. Svane, O. Gunnarsson, Phys. Rev. Lett. 65, 1148 (1990) 

9 Z.-X. Shen, R. S. List, D. S. Dessau, B. O. Wells, O. Jepsen, 
A. J. Arko, R. Barttlet, C. K. Shih, F. Parmigiani, Phys. 
Rev. B 44, 3604 (1991) 

10 V. I. Anisimov, I. V. Solovyev, M. A. Korotin, M. T. 
Czyzyk, G. A. Sawatzky, Phys. Rev. B 48, 16929 (1993) 

11 F. Manghi, C. Calandra and S. Ossicini, Phys. Rev. Lett. 
73, 3129 (1994) 

12 F. Aryasetiawan and O. Gunnarsson, Phys. Rev. Lett. 74, 
3221 (1995) 

13 M. Takahashi and J. Igarashi, Ann. Physik 5, 247 (1996) 

14 M. D. Towler, N. L. Allan, N. M. Harrison, V. R. Saunders, 
W. C. Mackrodt, E. Apra, Phys. Rev. B 50, 5041 (1994) 

15 B. H. Brandow, Adv. Phys. 26, 651 (1977) 

16 S. Hiifner, Adv. Phys. 43, 183 (1994) 



17 P. Fulde in Electron Correlations in Molecules and Solids, 
Springer Series in Solid State Sciences, vol. 100 (Springer, 
Berlin, 1993) 

18 J. Yamashita, S. Asano, J. Phys. Soc. Jpn. 52, 3514 (1983) 

19 H. Stoll, Phys. Rev. B 46, 6700 (1992); H. Stoll, Chem. 
Phys. Lett. 191, 548 (1992); B. Paulus, P. Fulde, H. Stoll, 
Phys. Rev. B 51, 10572 (1995); B. Paulus, P. Fulde, H. 
Stoll, Phys. Rev. B 54, 2556 (1996) 

20 K. Doll, M. Dolg, P. Fulde, H. Stoll, Phys. Rev. B 52, 4842 
(1995); K. Doll, M. Dolg, H. Stoll, Phys. Rev. B 54, 13529 
(1996) 

21 A. D. Becke, Phys. Rev. A 38, 3098 (1988) 

22 C. Lee, W. Yang and R. G. Parr, Phys. Rev. B 37, 785 
(1988) 

23 Note that, since localization is necessary, it will be more 
difficult to apply the method to metals, although an appli- 
cation to graphite has already been performed ( H. Stoll, 
J. Chem. Phys. 97, 8449 (1992)). 

24 R. J. Cave, E. R. Davidson, J. Chem. Phys. 89, 6798 (1988) 
W. Kutzelnigg in Methods of Electronic Structure Theory, 
edited by H. F. Schaefer III (Plenum, New York, 1977) 

26 R. J. Gdanitz, R. Ahlrichs, Chem. Phys. Lett. 143, 413 
(1988) 

27 H.-J. Werner and P.J. Knowles, Theor. Chim. Acta 78, 175 
(1990) 

28 see, e.g., G. D. Purvis III, R. J. Bartlett, J. Chem. Phys. 
76, 1910 (1982) 

29 K. Raghavachari, G. W. Trucks, J. A. Pople, M. Head- 
Gordon, Chem. Phys. Lett. 157, 479 (1989) 

30 MOLPRO is a package of ab-initio programs written by 
H.-J. Werner and P.J. Knowles, with contributions from J. 
Almlof, R. D. Amos, M. J. O. Deegan, F. Eckert, S. T. 
Elbert, C. Hampel, W. Meyer, A. Nicklass, K. Peterson, 
R. Pitzer, A. J. Stone, P. R. Taylor, M. E. Mura, P. Pulay, 
M. Schuetz, H. Stoll, T. Thorsteinsson, and D. L. Cooper. 

31 H.-J. Werner and P.J. Knowles, J. Chem. Phys. 82, 5053 
(1985); P.J. Knowles and H.-J. Werner, Chem. Phys. Lett. 
115, 259 (1985); H.-J. Werner and P.J. Knowles, J. Chem. 
Phys. 89, 5803 (1988); P.J. Knowles and H.-J. Werner, 
Chem. Phys. Lett. 145, 514 (1988); C. Hampel, K. Peter- 
son, and H.-J. Werner, Chem. Phys. Lett. 190, 1 (1992); 
P.J. Knowles, C. Hampel, and H.-J. Werner, J. Chem. 
Phys. 99, 5219 (1993) 

32 R. D. Cowan, D. C. Griffin, J. Opt. Soc. Am. 66, 1010 
(1976); J. H. Wood, A. M. Boring, Phys. Rev. B 18, 2701 
(1978) 

33 T. H. Dunning, Jr., J. Chem. Phys. 90, 1007 (1989); R. 
A. Kendall, T. H. Dunning, Jr., R. J. Harrison, J. Chem. 
Phys. 96, 6796 (1992) 

34 M. Dolg, U. Wedig, H. Stoll, H. Preuss, J. Chem. Phys. 
86, 866 (1987) 

35 CRC Handbook of Chemistry ad Physics, 75 th edition, Ed- 
itor: David R. Lide (CRC Press, Boca Raton, 1994/1995) 

36 H. B. Shore, J. H. Rose, E. Zaremba, Phys. Rev. B 15, 
2858 (1977) 

37 K. Schwarz, J. Phys. B 11, 1339 (1978), K. Schwarz, Chem. 
Phys. Lett. 57, 605 (1978) 

38 T. Grabo, E. K. U. Gross, Chem. Phys. Lett. 240, 141 
(1995) 

39 C. W. Bauschlicher and P. Maitre, Theor. Chim. Acta 90, 



5 



189 (1995) 

40 O. Madelung (editor), Landolt-Bornstein tables, Group III, 
Vol. 17g (Springer, Berlin, 1984) 

41 C. Wang, K. Fink, V. Staemmler, Chem. Phys. 192, 25 
(1995) 

42 F. Illas, J. Casanovas, M. A. Garcia Bach, R. Caballol, and 
O. Castcll, Phys. Rev. Lett. 71, 3549 (1993) 

43 C. de Graaf, R. Broer, and W. C. Nieuwpoort, preprint 
N. M. Harrison (Daresbury, UK), private communication 

45 R. L. Martin, P. J. Hay, J. Chem. Phys. 75, 4539 (1981); J. 
P. Desclaux, At. Data Nucl. Data Tables, 12, 311 (1973); 
P. Pyykko, Chem. Rev. 88, 563 (1988) 

46 R. A. Coy, C. W. Tompson, E. Gurmen, Solid State Com- 
mun. 18, 845 (1976) 

47 K.-H. Hellwege and A. M. Hellwege (eds.), Landolt- 
Bornstein tables, Group III, Vol. 7bl (Springer, Berlin, 
1975) 

48 S. Tanaka, J. Phys. Soc. Jpn. 62, 2112 (1993); J. Phys. 
Soc. Jpn. 64, 4270 (1995) 

49 S. Fahy, X. W. Wang, S. G. Louie, Phys. Rev. Lett. 61, 
1631 (1988); X.-P. Li, D. M. Ceperley, R. M. Martin, Phys. 
Rev. B 44 10929 (1991) 

50 C. E. Moore, Atomic Energy Levels, NSRDS-NBS 35 / Vol. 
I-III, Nat. Bur. Standards (Washington, DC) 

51 V. I. Srdanov, D. O. Harris, J. Chem. Phys. 89, 2748 (1988) 

52 L. R. Watson, T. L. Thiem, R. A. Dressier, R. H. Salter, 
E. Murad, J. Phys. Chem 97, 5577 (1993) 

53 D. W. Green, G. T. Reedy, J. G. Kay, J. Mol. Spectrosc. 
78, 257 (1979) 



FIG. 1. Spin ordering: a) in the ground state of NiO; b) 
for a double external configuration with simultaneous excita- 
tions from a Ni 2+ and an 2 ~ ion; c) for a single external 
configuration with spin-flip. 



FIG. 2. Result of the incremental expansion for NiO at a 
lattice constant of 4.17 A. 
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TABLE I. Parameters 

of the Ni 2+ pseudopotential. The pseudopotential has the 



form \f ( r 1 
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-0.2278919 


0.2415367 
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0.3984276 
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TABLE IV. Local correlation energy increments (in H) at 
a lattice constant of 4.17 A, from QDVPT calculations. The 
largest increments are compared with CCSD and CCSD(T). 
The numbers in brackets (a,b,c) are indices of the lattice site 
of the ion and are given in multiples of the lattice constant. 5 
denotes the spin coupling (5* = 2 means ferromagnetic, 5 = 
antiferromagnetic coupling) . The increments are already mul- 
tiplied with the weight factor except where the weight factor 
is explicitly repeated in the sum. A weight factor of means 
that this state does not occur in the AF2 state. In the case 
of the three-body increments, we summed only over the two 
largest ones in order to make a comparison with the results 
of Table V possible. 



TABLE II. Sum of first and second ionization potentials 
for Ni in Hartree units ( 3 F state) 





Ni -> Ni 2 + 


Hartree-Fock, finite differences 


0.840310 


Dirac-Fock, finite differences 


0.852950 


Hartree-Fock with relativistic pseudopotential, 




our basis 


0.855544 


QDVPT (3d, 4s correlated), our basis 


0.930863 


QDVPT (3s, 3p, 3d, 4s correlated), our basis 


0.936216 


LDA, our basis 


0.999031 


BLYP, our basis 


0.98256a, 


expt. 


0.94776CE] 



TABLE III. Bond length R e (A), dissociation energy 
D e (eV) and vibrational frequency ui e (cm -1 ) of the NiO 
molecule. The orbitals correlated with the multi-reference 
averaged coupled-pair functional (MR-ACPF) scheme are O 
2s, 2p, Ni 3d, 4s. The dissociation energy is calculated with 
respect to the 3 D (dV) state of Ni. The DFT calculations 
for the atoms were performed without applying a spherical 
approximation for the potential and with integer occupation 
numbers. 





Re 
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sum of one-body increments 


-0.026278 
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Ni(0,0,0)-O(2,2,l) 
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sum of Ni-O increments 


-0.073446 


Ni(0,0,l)-Ni(0,l,0), S=2 


3 


-0.000831 


Ni(0,0,l)-Ni(0,l,0), S=0 


3 


-0.000780 


Ni(l,0,0)-Ni(-1,0,0), 5=2 





-0.000025 x 


Ni(l,0,0)-Ni(-1,0,0), 5=0 
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-0.000025 x 3 


sum of Ni-Ni increments 


-0.001686 


O(0,0,0)-O(0,l,l) 
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-0.023496 


CCSD: -0.020136 






CCSD(T): -0.023928 
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O(0,0,0)-O(2,l,l) 
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-0.001200 


O(0,0,0)-O(2,2,0) 
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O(0,0,0)-O(3,l,0) 
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O(0,0,0)-O(2,2,2) 
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sum of O-O increments 
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O(0,0,0)-Ni(0,0,l)-Ni(0,0,-1),5=2 
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+0.002724 


O(l,0,0)-O(-l,0,0)-Ni(0,0,0) 
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+0.000078 


sum of two largest three-body increments 


+0.003572 


total sum 


-0.124189 


expt. 


-0.147 
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TABLE V. Local increments (in H) for NiO at a lattice 
constant of 4.264 A. 





QDVPT | 


sum of one-body increments 


-0.028014 


sum of Ni-O increments 


-0.068362 


sum of Ni-Ni increments 


-0.001416 


sum of O-O increments 
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sum of three-body increments 
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total sum 
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